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Abstract 


In this paper, we introduce a modified interpretation of the water activity presented in Springer et al. [T.E. Springer, T.A. Zawodzinski, S. 
Gottesfeld, Polymer electrolyte fuel cell model, J. Electrochem. Soc. 138 (8) (1991) 2334-2342]. The modification directly affects the membrane 
water transport between the anode and the cathode (two electrodes) of the polymer electrolyte membrane (PEM) fuel cell in the presence of 
liquid water inside the stack. The modification permits calibration of a zero-dimensional isothermal model to predict the flooding and drying 
conditions in the two electrodes observed at various current levels [D. Spernjak, S. Advani, A.K. Prasad, Experimental investigation of liquid 
water formation and transport in a transparent single-serpentine PEM fuel cell, in: Proceedings of the Fourth International Conference on Fuel 
Cell Science, Engineering and Technology (FUELCELL2006-97271), June 2006]. Using this model the equilibria of the lumped water mass in 
the two electrodes are analyzed at various flow conditions of the stack to determine stable and unstable (liquid water growth) operating conditions. 
Two case studies of water management through modification of cathode inlet humidification and anode water removal are then evaluated using this 
model. The desired anode water removal and the desired cathode inlet humidification are specified based upon (i) the water balance requirements, 


(ii) the desired conditions in the electrodes, and (iii) the maximum membrane transport at those conditions. 


© 2006 Published by Elsevier B.V. 
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1. Introduction 


Due to the significant influence of water [3,4] on the perfor- 
mance of the polymer electrolyte membrane fuel cell (PEMFC), 
many models have been developed to capture its spatial distribu- 
tion within a fuel cell stack. These computational fluid dynam- 
ics (CFD) models are then used to optimize the cell flow field 
(width, length, orientation, etc.) and materials (hydrophobicity, 
porosity, thickness, etc.). Once the optimum fuel cell is designed, 
CFD models [5] can again, theoretically, be used to select the 
right number of cells to be combined in a stack in order to avoid 
large cell-to-cell variations. Although the above fuel cell design 
process is visionary, it is still far from being realized. Ultra large- 
scale models [6] have so far been used to simulate medium size 
single cells with fixed geometry and materials interconnection. 
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Once the fuel cell is designed the only remaining degrees of 
freedom for water management are the mass flow rates entering 
and leaving the electrodes of the fuel cell stack. The actuating 
mechanisms for controlling water flow rates, for example, anode 
and cathode inlet humidification, oxygen excess ratio (OER), 
and anode recirculation, are not spatially distributed, but allow 
control only in the stack manifolding, which in turn constitutes 
the boundary control to the electrodes. Water management of 
the fuel cell stack with these actuators cannot be used to control 
the spatial distribution of water in the stack (from cell to cell 
and within the stack material). Nevertheless, one can start by 
controlling the mass of water in the two lumped volumes that 
correspond to the anode and the cathode electrodes and aug- 
ment this methodology with diagnostic mechanisms that detect 
maldistribution of water through cell to cell variability or large 
pressure variations [7]. It is through the water management in 
the electrodes that we hope to control the humidity conditions 
at the catalyst where the reaction happens. 

In this paper, the capabilities of physically actuated system 
variables to achieve water balance in the stack are evaluated 
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Nomenclature 

p pressure (Pa) 

m mass (kg) 

y vapor mass fraction 
W flow rate (kg s7!) 

T temperature (K) 

Igt stack Current (A) 


N molar flux (mole s7! cm7?) 
M molecular weight (kg mole~!) 
R gas constant (Jkg~! K7!) 

À water content 


dw water activity 

(0) relative humidity 

w humidity ratio 

k flow coefficient (kg s7! Pa!) 


Superscripts and subscripts 


a anode 

c cathode 

m membrane 
st stack 

v vapor phase 
l liquid phase 
H20 water 

H2 hydrogen 
O2 oxygen 

N2 nitrogen 
gen generated 
rem water removed 
in inlet 

out outlet 

sat saturated 


using a zero-dimensional, isothermal, lumped parameter 
PEMFC stack model similar to [8]. In this low order model, wa- 
ter (both liquid and vapor) accumulates in the lumped volumes 
of the two electrodes. Additional components in the stack that 
can potentially store water, such as gas diffusion layer (GDL) 
and membrane, are neglected. Moreover, the model assumes in- 
stantaneous transition from vapor to liquid phase if the partial 
pressure of water vapor exceeds the saturation pressure. Finally, 
an ideal thermal control is assumed that maintains constant tem- 
perature throughout the stack. 

Clearly, the stack model with these assumptions will not cap- 
ture any spatial distribution. Instead, this low order model fo- 
cuses on important lumped phenomena, such as, anode flooding 
at lower current densities [2,9] and anode drying at higher cur- 
rent densities [10]. In order to capture the phenomenon of anode 
flooding, the model presented in [12,11] is modified to allow 
liquid water on the cathode side of the membrane to influence 
membrane water transport. 

The low order model facilitates the stability analysis of the 
water equilibria arising from employing various popular water 
management techniques, such as controlling humidity of inlet 


flow or anode water removal using recirculation [14,13]. This 
step is important since only stable equilibria can be achieved 
via feedforward maps (look-up tables) that use the information 
of the current drawn to determine, for instance, the appropriate 
level of anode water removal [15] or cathode inlet humidification 
[16,17]. 

We demonstrate a methodology to achieve water manage- 
ment using (i) cathode inlet humidification, and (ii) anode water 
removal, as actuating mechanisms. Other mechanisms, namely, 
anode humidification and cathode water removal through high 
oxygen excess ratio, might be necessary when operating at very 
low or very high currents. The use of these mechanisms will 
need more investigation than the one that can be allotted in this 
paper. 

The paper is structured as follows: the stack model related 
to water management is discussed in Section 2. The water 
dynamics in the anode and the cathode and their interaction 
under different operating conditions are discussed in Section 
3. The model predictions of stack flooding and drying char- 
acteristics without active water management are discussed in 
Section 4. In Section 5, different water management strate- 
gies involving cathode inlet humidification and anode water 
removal are identified and their ability to achieve particular 
water conditions are investigated independently of specific ac- 
tuating components. Section 6 summarizes the results in the 


paper. 
2. Model for water management analysis 


In this section, we first present the membrane water transport 
model [1] with a modification to the typical interpretation of 
water activity used in [11,12,18]. The proposed modification 
enables the prediction of liquid water accumulation in the anode 
at low current densities, which is important for mobile fuel cell 
applications due to the frequent operation at low loads. Then, 
we integrate this modified membrane water transport model with 
anode and cathode mass balance equations to develop the two- 
volume fuel cell model. 


2.1. Membrane water transport 


The water transport through the membrane from the anode 
to the cathode, Wmemb, is given by 


Wmemb = Myo Afelcetts(Narazg — Naite), (1) 


where Age(cm?) denotes the active fuel cell area, nels the num- 
ber of cells in the stack and Mmo the molecular weight of water. 
The molar flux per unit cell from the anode to the cathode due to 
the electro-osmotic drag is Narag (mole s7! cm~’), while Naitf 
(mole s7! cm~?) is the molar flux per unit cell from the cathode 
to the anode due to the difference in water concentration on the 
anode and cathode sides. 

These fluxes depend upon the stack current, Ist; the membrane 
water content, Àm; and the amount of water present on the anode 
and cathode sides [1,11]. The molar flux due to electro-osmotic 
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drag from anode to cathode, Narag, is given by [1] 


N =j Ist (2) 
drag = d An F , 


with 
na = 0.002942, + 0.05Am — 3.4 x 107}, 
where F is Faraday’s constant. 


The molar flux due to back diffusion from cathode to anode, 
Naite, is modeled as [1] 


Pm,dry 
Naitt = «p Dw ——— 


(Ac — Aa), (3) 
tm Mm,dry i i 


with 


Io 1 
D, =D xal oo 
WT A EAR ( (= ~)) 


where Tst is the stack temperature, tm is the membrane thickness, 
Pm,dry and Mm,ary are the density and the weight per mole, re- 
spectively, of the membrane when dry. The diffusion coefficient, 
Dy, is scaled using a multiplicative term ap to serve as a tuning 
parameter. In our analysis we use œp in evaluating the effect of 
membranes with different diffusion properties.! The coefficient, 
D,, is defined in [12,11] as 


1076 if Am <2 
1076(1 +2(àÀm —2)) if2<ìm<3 
1076(3 — 1.67(Am — 3)) if 3 < Am < 4.5 
1.25 x 1076 if Àm > 4.5. 


D (4) 


In the above equations, A, and A, are the water contents at the 
cathode and anode side of the membrane, respectively. Water 
content [1,19] is defined as 


Ay) = flaw); (5) 


where dw .) represents the water activity and the function, f} is 
shown in Fig. 1. Here, subscript “(-)” is to be replaced with “a”, 
“c” and “m” for anode, cathode and membrane, respectively. 

Two important characteristics of the relation between dy, (.) 
and A,.) are that it is monotonic and has an upper bound. The up- 
per bound on the membrane water content, Amax, corresponds to 
the water content of the membrane when completely immersed 
in liquid water [1]. In this analysis, this maximum water content 
is achieved once the water activity equals 3, beyond which the 
water content is constant [19]. 

The average membrane water activity, dw.m, iS defined? as 


awa + aw.c 
a ae, 
where awa and awc correspond to the water activities on anode 
and cathode sides of the membrane, respectively. 


(6) 


adw.m = 


1 The term ap can also be combined with D}, thereby making D; as the tuning 
parameter. 

2 Note that in the following work aw,m is defined as the average membrane 
water activity, similar to [12]. Other possible definitions that could be employed 
[20] are aw,m = aw,a to reflect the importance of the anode side water activity, 
Or dw,m = Min(dw,a, Gw,c) to capture the conditions of the limiting electrode. 


18 T 7 T T r 7 
16 4 =16:8 
147 Linear approximation 7 
Springer et al. [1] Water activity 
= 12) Subsaturated does not affect | 
© „| conditions water content 
© 10 : 
= 
fo] 
O g 4 
o 
ta Saturated conditions 
i 7 
4 H oa 
2 al 
% 0.5 1 1.5 2 2.5 3 3.5 4 


Water activity (a,, ) 


Fig. 1. Water content as a function of water activity. 


In this simplified isothermal model, the water activity on the 
anode or the cathode sides of the membrane is defined as 
aw.) = ae (7) 
17H20.v.(-) 
where my,0,(.) is the total water in the adjacent electrode (anode 
or cathode). The maximum water vapor that can be present in 
that electrode, miso vc) 1S given by 


sat, Psat (To) Yo 


™yH,0,v,(-) = Ry,oTo (8) 


where Psat(7(.)) denotes the saturation pressure corresponding to 
temperature T), where Ti) is the temperature, Vq. is the lumped 
volume of the electrode and Ry,0 is the gas constant of water 
vapor. 

With ay, (.) defined as Eq. (7), Wmemb is influenced by liquid 
water in the stack and the isothermal stack model equations used 
here can capture anode flooding behavior under conditions ex- 
perimentally observed [21,2], when the diffusion coefficient, ap, 
is appropriately tuned or defined. Different modifications can be 
found in literature [19,22] that incorporate effect of liquid water 
on membrane water transport, however, these are applicable to 
more detailed models. 

Note here that, modeling adw,.) to be equal to the relative 
humidity in the adjacent electrode, ¢,.), as in [18,11,12], instead 
of Eq. (7), always results in sub-saturated anode conditions even 
when the cathode is flooded. Indeed, if 


PH20,v,(-) 


f 9 
Psat(T(.)) A 


aw) = PO) = 
where pp,o,v,() is the partial pressure of water vapor in an 
electrode, then the water activity is always less than 1, as 
PH30,v,(.) < Psat(7(.)). Even when liquid water is present on the 
cathode side, using the activity definition (9), we have ay. = 1. 
With the maximum water activity in the cathode being 1, the 
water transported from cathode to anode is limited and can at 
most cause saturating condition in the anode. Specifically, when 
the partial pressure of water vapor in anode reaches saturation, 
awa = 1, the water transport from cathode to anode stops be- 
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(Pin,as din,a) (Pinc Pin,c) 


(Pout,a) 


(Dout,c) 


Fig. 2. Lumped volumes representing the two electrodes of the fuel cell stack. 


cause the diffusion requires a positive gradient in water content 
between anode and cathode, that is, Ag < Ac. This implies that 
an isothermal stack model using water activity definition (9) will 
never exhibit anode flooding. 

The model for water mass transport through the membrane 
is next used in the mass conservation equations applied to the 
lumped volumes of the anode and cathode electrodes of the fuel 
cell. The resulting ordinary differential equations are then used 
to find the mass equilibria and stability properties for the dry gas 
and the water inside the two electrodes. 


2.2. Mass conservation 


The fuel cell stack equilibria are derived by assuming isother- 
mal conditions in the electrodes and using mass conservation for 
the important species shown in Fig. 2, namely, the water in the 
anode and the cathode, my,0,q and my,0,c, respectively, the hy- 
drogen in the anode, my, a, and dry air in the cathode, mo,Nn,,c- 
The water mass in the two electrodes are given by [12] 


dm, 0.2 
2U,a 
a = Wu,0,in,a Sane Wu,0,out,a = Wmemb (10) 
dmp o 
2U,C 
—ap = We 0,gen F Wu,0,in,c = Wu,0,out,c + Wmemb, 


(11) 


where W4,0,in,a and Wy,0,out,a are the mass flow rates of water 
entering and leaving the anode, respectively. Similarly, W4,0, in,c 
and Wy,0,out,c are the mass flow rates of water entering and 
leaving the cathode. The water generation rate due to the elec- 
trochemical reaction is given by 

I st 


WH30,gen = Neells 7 p Mto- (12) 


The mass flows in the two electrodes can be derived through 


generic equations, where “(-)” should be replaced by “a” for an- 
ode and “c” for cathode. The mass of water entering an electrode 


is given by 
Web, in,-) = Yv,in,(-) Win, (13) 
Win.) = Kin) (Pin) — PO), (14) 


where yy,in,.) is the mass fraction of water vapor in the inlet 
flow. The total flow entering an electrode, Win,.), is modeled 
using a linear flow model with a flow coefficient kin, Here, 
Pin,(.) is the total pressure at the electrode inlet while pv.) is the 
total pressure within the electrode. 

The mass fraction of water vapor in the inlet flow is 


PH20,y, in.) MH, 0 


Yv,in,(.) = T 7 i 
- PH0,v,in,.) Mmo + (Pin) — PH20.v,in)) Mg, 
(15) 
with the partial pressure of vapor at the electrode inlet as 
PH0,v,in,) = Qin, Psat(Tin,C)), (16) 


where Psat(Tin,(.)) denotes the saturation pressure corresponding 
to inlet temperature Tin,(), and ¢in,(.) is the relative humidity of 
the inlet flow. Here, we assume that no liquid water is supplied 
to the electrodes. 

Meanwhile, M is used to represent molecular weight and 
the subscript “g, (-)” denotes the dry gas reactant in the elec- 
trode, and hence, is replaced by “H2” for hydrogen in anode and 
“O.N>” for dry air in the cathode. 

Using Dalton’s law of partial pressures, the total pressure in 
an electrode is given by 


PC) = PHO,v, 0) F Pg); (17) 


where pg.) and pp,o,v,¢) are the partial pressures of the dry gas 
and the water vapor in the electrode. 

Considering the ideal gas law, py,o0,v,(.) depends upon the 
mass of water vapor in the electrode, my,0,y,(.), a8 


= —>— mH, 0,()- (18) 


Assuming instantaneous condensation and evaporation of water, 
the mass of water vapor is given by 


M¥H0,v,(-) = min {MH30,(.); mitov o? , (19) 


sat ; ; 
where jy, 0, y,.) 18 the maximum water vapor that can be present 


in the electrode (see Eq. (8)). 

The partial pressure of dry gas, pg,.), is obtained from the 
mass of the dry gas in the electrode, mg_.), using the ideal gas 
law, as 
Pa) = =I g,(.); (20) 
where Rg (.) represents the gas constant for the dry gas in the 
electrode. 

The water leaving an electrode is separated into liquid and 
vapor phases due to differences in their flow properties. 


Wep0, out,(-) = WH» 0,v,0ut,(.) + WH20,1,0ut,()> (21) 


where W4,0,v,out,(-) and WH, 0,1 out,(.) are the flow rates of water 
vapor and liquid, respectively, leaving the electrode. 
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The flow rate of water vapor flowing out of an electrode is 
given by 


We 0,v,out,(-) = Yv,(-) Wout, (22) 
Wout,(.) = kout OPO — Pout,()), (23) 


where Pout,(.) is the pressure at the electrode outlet, Wout.) is the 
net mass flow rate of the gases leaving the electrode, and Kout,(-) 
is the linear flow coefficient for outlet flow. The mass fraction 
of water vapor in an electrode is 
™¥H0,v,(-) 
My0,v,() + Mg.) 


Wi) = (24) 
Based on [23], the flow rate of liquid water drained out of an 

electrode, Wy,0,1,out,(.), depends upon the accumulated mass of 

liquid water and the pressure drop across the electrode as 


Wu0,1out,.) = (PO — Pou, ©) f 0Nnm 0,1,0), (25) 


where function f(mp,o,1,()) depends upon channel geometry 
and material properties. The mass of liquid water accumulated 
in the electrode, my,0,1,.), 18 given by 


MH0,1,(-) = MH0,() — 1H 0,v,(-)- (26) 


In order to complete the stack model, we need to provide 
the mass of dry gas in the anode and cathode electrodes. The 
dynamics of mass of hydrogen in the anode, my, a, are modeled 
by applying conservation of mass, as 


dmy, a 
dt 


where Wy, ina and Wp, out,a are the flow rates of hydrogen at 
anode inlet and outlet, respectively. The rate at which hydrogen 
is consumed in the electrochemical reaction is given by 

I st 


Why, ret = Necells zp Mit: (28) 


Similarly, the dynamics of mass of dry air in the cathode are 
given by 


= Why, in,a = Wh;,out,a = We, ret, (27) 


dmonNo,c 
dt 
where Wo Nj, in,c and Wo,N3,out,c are the flow rates of dry air at 
cathode inlet and exit, respectively, and the rate at which oxygen 
is consumed in the electrochemical reaction is given by 
Wop, ret = edie Mes (30) 
i 4F 
where Mo, is the molecular weight of oxygen. 
The flow rates of dry gas into and out of the electrode are 
given by 


We in) = (l — Win) Win, © (31) 
Wg, out,(-) =(1- Yv,¢-)) Wout, > (32) 


where Yy,in,() and yy out.) are the vapor mass fractions in the 
inlet and exit flows, Win.) and Wout,.), respectively. 

Oxygen and nitrogen in the cathode are not modeled sep- 
arately, since we are interested only in the total pressure and 
water dynamics in the electrode. These dynamics depend upon 


= WowN;, ine — WO Nd, out,e — Wo, ret (29) 


the molecular weight of the gas mixture which is assumed con- 
stant at29 x 107°? kg mole™!. Even though oxygen is consumed 
inside the cathode, this assumption introduces, at the most, a 4% 
modeling error in the molecular weight of the gas mixture, as 
the molecular weight of the gas mixture is bounded between 
28 x 1077 and 29 x 10-7 kg mole™!. 

Eqs. (10), (11), (27) and (29) obtained through mass conser- 
vation are used by the two-volume model to predict the mass 
equilibria in the lumped volumes of the two electrodes. In these 
equations, the electrode inlet conditions, that is, inlet pressures 
Pina and pin,c, and inlet humidities ģin,a and din, modify the 
flows into the two electrodes. Meanwhile, the flows out of the 
electrodes are affected by the electrode exit conditions Dout,a 
and Pout,c- The different inputs and dynamic states of the two- 
volume model can be summarized using Fig. 2. 

We show next a stability analysis of different water mass 
equilibria predicted by this simple two-volume model (Section 
3). In Section 4, we state the model predictions of flooding 
and drying conditions along with a discussion of its key cal- 
ibrating parameters and we show how the model predications 
can be used to analyze the suitability of water management 
schemes. 


3. Analysis of equilibrium conditions related to water 
management 


This simplified model can describe complex phenomena as- 
sociated with the water mass in the stack. For example, consider 
the case when the inlet and exit conditions of the stack are cho- 
sen such that the system is at an equilibrium. If this equilibrium 
condition is unstable, then in the absence of a proper water man- 
agement scheme a small change in the stack current can lead to 
an unbounded growth in the mass of water inside the stack, as 
shown in Fig. 3. 

As time progresses the mass of water in the cathode reaches 
a level such that the cathode water activity exceeds 3, seen in 
Fig. 3 at 10s. Once this condition is reached, the additional 
water accumulating in the cathode does not increase the back 
diffusion from cathode to anode, thereby allowing the mass of 
anode water to equilibrate at a value depending upon the anode 
inlet and exit flows and the electro-osmotic drag. Thus, Fig. 3 
reflects the different stability properties of the anode and cathode 
water, in presence of a disturbance from stack current. 

In this section, equilibrium conditions at which the mass of 
water in the anode or cathode can exhibit an unstable behavior 
are identified. First, physical reasoning is used to justify the 
dynamic behaviors at different operating conditions. Then, a 
linearized approximation of system dynamics is used to verify 
the results analytically. 


3.1. Stability of different operating conditions 


Stability of typical operating points corresponding to differ- 
ent conditions are listed in Table 1. The columns (a) and (b) 
address the effect of liquid water removal on the stability of the 
water dynamics in the system. In each case, the behavior of mass 
of water in anode and cathode is specified separately. 
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Table 1 
Stability of the water dynamics under different operating conditions 


Case Water level Stable equilibria* 
Be ; aw dW, 
Ba Zoa (a) H20, 1, out, (-) =0 (b) H20, |, out, (-) +0 
Mi,0,v,c Mi O.v,a dMH0,1,.0) dMH0,1,0) 
Cathode Anode Cathode Anode 
(i) [0 1] [0 awc] Yes Yes Yes Yes 
(ii) [1 3] [0 1] Yes Yes Yes Yes 
(iii) [1 3] [1 awe] No No Yes Yes 
(iv) >3 [0 1] No Yes Yes Yes 
(v) >3 [13] No Yes Yes Yes 
Siinode the pressure, p;.), and mass fraction of water vapor, yy ,(.), inside 


Time (sec) 


Current (A) 


0 5 10 15 
Time (sec) 


70 = 


Fig. 3. Unstable water dynamics in the stack without water management: a 
response to a change in current drawn. 


The rows in Table 1 distinguish the operating condi- 
tions in the stack based upon the water levels in the an- 
ode and the cathode electrodes. We first differentiate between 
sub-saturated (my,0,(.)/ mito, o< 1) and saturated conditions 


(mn0,0/ MiTo, > 1) due to the dissimilar effects of the wa- 
ter vapor and liquid water on the electrode pressures. Within 
the saturated conditions, we also differentiate between satiated 
d< MH 0,()/ Myo.) < 3) and soaked (m¥,0,()/™ 50.) > 
3) conditions due to differences in their water content re- 
sponse (see Fig. 1), and hence the water diffusion through the 
membrane.* Due to space limitation, Table 1 shows only the 
conditions for which aw, a < aw,c, and a few of the conditions 
listed are elaborated here. 

Consider case (i) where the water in both the anode and the 
cathode is only in vapor phase and a water balance is achieved 
in the anode and cathode. The interactions between the differ- 
ent states, that is, the lumped mass of cathode water, my,0.c, 
anode water, my,0,a, anode hydrogen, my,,q, and cathode dry 
aif, M0,N>,c, are shown in Fig. 4. The dash-dot arrows indicate 
the flows in and out of the electrodes and the membrane water 
transport. The flows in and out of the electrodes are affected by 


3 More sub-categories can be established by considering the various possible 
values or models for the membrane water content. We will not elaborate on the 
anomalies that depend upon the definition membrane water activity. 


the electrodes, as seen in Eqs. (13), (21), (31) and (32). Both 
Pv.) and yy (.) depend upon the mass of water vapor and the mass 
of dry gas in the corresponding electrode, and this interaction 
is shown by the solid arrows in the figure. It is the coupling be- 
tween the mass of a species in an electrode to the different inlet 
and exit flows that defines the stability of the equilibrium. 

For example, consider a change in the stack current from its 
equilibrium value for a short duration that generates additional 
water in the cathode and disturbs the water balance. Although the 
stack current returns to its original value, the increase in cathode 
water leads to an increase in the cathode humidity. Consequently, 
the mass of water in the anode increases due to an increase in 
diffusion through the membrane. The increase in anode humidity 
will reduce the gradient across the membrane, and thus reduce 
the diffusion. Meanwhile, the increase in the anode and cathode 
humidities will increase the exit water flows, thereby causing 
the mass of water in the anode and the cathode to return to their 
original levels such that the exit flows balance the inlet flows. 
This self-balancing mechanism between the electrode humidity 
and the water transport across the membrane or in/out water 
mass flows establishes the stability of the operating points in 
case (i). 

Now consider a case where the dynamics of both anode and 
cathode are unstable, for example, case (iii)-(a), where the liq- 
uid water does not leave the stack through the exit flows even 


Memb. water transport 
` 


774 (Aa) 


1 we-8 


Water flow 
(in/out) 


Water flow 
(in/out) 


adw aiis 
la a 
1 3 oQ ! 
1 = = il 
I à S 1 
i “EN ea | 
1 gl Q I 
1 I Qi 
` Di tt 
“1° “T° 
I I 
N essa 
Hydrogen Dry air 
flow (in/out) flow (in/out) 


Fig. 4. Interactions between different states in the system when anode and cath- 
ode do not contain liquid water. 
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Fig. 5. Interactions between different states in the system when both anode and 
cathode contain liquid water and Wy, 0,1,out,(.) = 0. 


though liquid water is present in the anode and the cathode. The 
interaction between the various states of the system in this case 
is shown in Fig. 5. Both the pressure and the mass fraction of 
water vapor depend upon my,0,y,(.), Which remains constant at 
its maximum value Mi Ov,()" The broken lines shown in the 
figure imply that the mass of water in an electrode does not af- 
fect the inlet or exit flows due to their inability to affect pọ and 
WC): 

In this case, a perturbation in the mass of water inside the 
anode or cathode does not affect pi.) or yy,.), and hence, there is 
no change in the flows entering or leaving the stack. Note that, 
the total mass of water inside the stack cannot be affected by the 
membrane water transport. There is no stabilizing mechanism 
to restore the states back to their original value once they are 
perturbed. The mass of water in the anode and cathode electrodes 
may move to a different operating point depending upon the 
perturbation. 

There are cases, such as the one shown in Fig. 3, that even 
though the mass of water in the cathode increases continuously 
due to the additional water being generated, the increase in an- 
ode water level ceases after some time (10 s in the case of Fig. 3). 
This phenomenon can be explained by analyzing the dynamics 
of an equilibrium condition in case (v)-(a), where the mem- 
brane water transport is independent of the mass of water in the 
cathode. 

Once the cathode water content attains its maximum value, 
Ac = Amax (see Fig. 1), the diffusion of water from cathode to 
anode becomes independent of the mass of water in the cathode. 
In Fig. 6, the broken line between my,0,- and the membrane 
water transport indicates that the membrane water transport does 
not change due to a perturbation in my,0,c. The interactions 
between the different states of the system indicate that my,0,c 
does not affect the mass of water in anode. The membrane water 
transport depends upon my,0,a, and contributes to the stability 
of mMH,0,a. 

The various stability results explained here using physical 
reasoning were validated using simulations. The results, how- 
ever, can be verified analytically by using linear approximation 
of the system dynamics under small perturbation. 
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Fig. 6. Interactions between different states in the system when anode and cath- 
ode contain liquid water and the membrane water transport is independent of 
mass of water in cathode. 


3.2. Analysis using linearization of system dynamics 


Consider the fuel cell system at equilibrium with the mass 
of water on the cathode and anode sides as MiH,0, < and mi,0, m 
respectively, the mass of hydrogen in the anode as mh, as and 
dry air in the cathode as M6,N>,c* Assume that the states are 
perturbed from their equilibrium values by x4,0,c, XH20,a, XH3,a 
and xo,N3,c, respectively, such that 


Mi.) = MF) + Xi), (33) 


where “i, (-)” is to be replaced with “H20, c”, “H20, a”, “Ho, a 
and “O2N2, c” for water in the cathode, water in the anode, 
hydrogen in the anode, and dry air in the cathode, respectively. 

The linear approximation of the dynamics of the perturbed 
system with states XH,0,c» XH20,a> XHy,a aNd XO,N),c, is obtained 
by using a Taylor series expansion of Eqs. (10), (11), (27) and 
(29) and neglecting the higher order terms. The linearized dy- 
namics of the perturbed system are represented using a state- 
space matrix notation [24] as 


XH)0,c XH20,c 
d | xH0,a ae, XH)0,a (34) 
dt XH3,a XH3,a 
XO2N2,c XO2N2,c 
with 
= ky0,c Ma 0 KO5Np,c 
Me [Me —Ma — kH,O,a kb, a 0 
0 —4H0,a qUp,a 0 
—qH20,c 0 0 —4ONd,c 
(35) 


where u(), ki) and qi.) are the coefficients derived from the 
linearization of the membrane water transport, flow of water 
into and out of the electrode, and gas flow into and out of the 
electrode, respectively. These coefficients depend upon the op- 
erating conditions associated with the equilibrium about which 
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the system is perturbed, and are important for establishing the 
stability properties of the perturbed system (34). 

The coefficients related to the flows of water in the electrodes 
are given by 


OW, 0,1, out,(-) 


OW4H,0,v, out, (-) 


OW4p0, in, (-) 
Omi.) 


ki) = 


(36) 


and after substituting WH, 0, in,.), WH20,v,out,(.) and WH 0,1,out,(-) 
based on Eqs. (13), (22) and (25), respectively, and taking into 
account that the inlet pressures, pin .), and water vapor mass 
fraction, yy,in,(.), depend solely upon electrode inlet conditions 
and not on the electrode states, we get 


dpc.) 
dmi) 


dpc.) 


ki.) = —Yv,in,()Kin,() (- ) + Yv, O kout, ©) 


ƏYy,( OW 0,1, out, 

+ Wout.) 2 a 
Omi.) Əmi) 

: 3 OPC.) 
= Wv,in,(-)Kin,(-) + Yv, k t,(- a 

(Yv,in, Okin, © v,() out,(.)) ðmi 
ƏYy,( OW H7O,1,out,( 

+ Wout, () WO os . (37) 
Əmi) Əmi, 


Similarly, the coefficients corresponding to the dry gas flows 
in the electrodes are given by 


OWe inc AW, l 
qio = -EmO + SEY, (38) 
Əmi) Əmi) 
and after using Eqs. (31) and (32) in Eq. (38) we get 
7 AP) 
to = -0 = wining (- 
i,(-) 'v,in,(-))*in, (-) ami. 
i APC) Iv.) 
1 — yy, © Kout, : Wout.- 
+( Wv,( )) out, (-) miO ae ami.) out, (-) 
5 z Op. 
= ((1 — yin, kine) + CL — yv )kou, o) 
3Y.) 
—W. j——. 39 
out, (-) nee (39) 


Eqs. (37) and (39) show that the coefficients kj.) and qj,.) 
depend on how the pressures and water vapor mass fractions in 
the electrode, and the liquid water removed from the electrode, 
are affected by the accumulation of mass of species “i, (-)”, that 
is, dp.) /0Mj,(.), OYy,¢.)/OMj,(.) and OWy,0,1,out,(.)/OMi,(-), respec- 
tively. 

The partial derivatives of pressure and mass fraction of wa- 
ter vapor in an electrode with respect to the dry gases in that 
electrode are given by 


Po) _ Re oTo 


as (40) 
OM g .(.) Vo 
ƏYv, e) = Yv,(-) (41) 
Mg.) MH0,() + Mg,(.) 


where “g, (-)” is replaced by “Ho, a” for hydrogen in anode and 
by “O2N2, c” for dry air in the cathode. 


The partial derivatives of pressure and mass fraction of water 
vapor in an electrode with respect to water in that electrode are 
given by 


QP) _ RmoTo 3mm o.v.) (42) 
dM H0.(-) Vo 3mm o,o) 

awo LW. mH O.v,0) (43) 
dmm oO MMO, + Mg) MMO) ` 
where 
ðmmov o _ J 1 ifmmo,v,© < Mihov.) (44) 
IMH)0,(.) 0 otherwise. 


The equation above indicates that we need to evaluate the sta- 
bility of the system separately for saturated and sub-saturated 
conditions inside an electrode. 

Finally, the coefficients related to membrane water transport 
are given by 


_ aW memb = OWmemb 


Ha A 
ƏMH,O,c 


ƏMHO,a 
where Wmemb is given in Eq. (1). Note that, u) depends upon 
the relation between lumped mass of water in the electrode and 
the corresponding water content, A,.) (see Fig. 1). 

The linearized model, thus obtained, is used to determine 
the stability characteristics of an equilibrium condition [24], by 
calculating the eigenvalues of the matrix A* given in Eq. (35). 
The eigenvalues, 5, are given by the solution to the characteristic 
equation 


det{514,.4 — A*} = 0, (45) 


where 74x4 is an identity matrix. The system is indeed stable 
when all the four eigenvalues, 5, have negative real parts. 

Due to the sparse nature of matrix A*, it is possible to re- 
arrange the rows and columns and isolate the states associated 
with the stable and unstable eigenvalues. For example, in case 
(iii)-(a), with liquid water present in both the electrodes, A* 
given in Eq. (35) can be simplified to 


—He Ha 0 —koyno,c 
— — kp, 0 
At Uc Ha Hp,a (46) 
0 0 —qH;,a 0 
0 0 0 —ON>,c 


This simplification is obtained by noting that dp(.)/dmy,0,.) = 
0, dYy.)/OMH,0,(.) = 0, and OW4,0,1,out,(.)/077H30,(-) = 0, 
which leads to 


kmo) =0, qmo, =0. 


In this case, A* has an upper “block” triangular form, hence, the 
dynamics of the dry gases xH, ,a and x0,N,,c are independent of 
the water dynamics and exhibit stable behavior, since gy, ,a and 
O.N>,c can be shown to be strictly positive by substituting the 
derivatives given by Eqs. (40) and (41) in Eq. (39). 

The remaining coefficients in .A* correspond to the water dy- 
namics in the two electrodes. The linear dependency between 
the first and second columns indicates that one of the eigenval- 
ues is zero and the subsystem stability can only be determined 
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through the analysis of the nonlinear system dynamics.* Indeed, 
the nonlinear system described with the total mass of water in 
the two volumes by adding Eqs. (10) and (11) is governed by 
integrator-like dynamics. Therefore, once perturbed, the mass 
of water in the anode and the cathode will deviate from their 
equilibrium conditions in case (iii)-(a), leading to the instability 
conclusion. 

For the case (v)-(a), the coefficient pe = 0, since the cathode 
water activity exceeds 3, and hence, the membrane water trans- 
port now depends only upon mass of water in the anode, which 
leads to 


0 Ha 0 —kooNn>,c 

At = O Ha —kyy,a 0 (47) 
0 0 —qH;,a 0 
0o 0 0 —4ONp,c 


with the eigenvalues of — ua, —GH),a, —GO,N>,a and 0. As dis- 
cussed for case (iii)-(a), the eigenvalue at the origin is due to the 
integrator corresponding to the total mass of water in the stack, 
and is not an artifact of linearization. Meanwhile, the eigen- 
value — ua, with ua > 0, corresponds to xy,0,a, indicating that 
the mass of water in the anode is stable even when the total mass 
in the stack grows unbounded. 

The discussion of these few cases is not exhaustive but serves 
in highlighting the complexity of the water dynamic behavior 
that the simple two-volume model can predict. The predictions 
of this two-volume model are now presented under certain in- 
teresting conditions. First, the predictions without water man- 
agement are discussed, and are then followed by the predictions 
with water management using cathode inlet humidification and 
anode water removal. 


4. Model simulation with no water management 


In this section, the steady state stack hydration characteristics 
without appropriate water management are predicted using the 
two-volume model for various fuel cell loads. The influence 
of the membrane diffusion parameter on these patterns is also 
considered in order to illustrate how the model can be tuned 
to match experimentally observed patterns. The prediction of 
flooding and drying patterns is useful in identifying the water 
management strategies required at different loads. 

The predictions of the two-volume model depends on the 
choice of the inlet and exit pressures. Physically realizable inlet 
and exit conditions are required to predict meaningful hydration 
characteristics at different fuel cell loads. In this paper, the two- 
volume model is simulated using the stack operating conditions 
of a high pressure fuel cell system. 

The important features of the operating conditions chosen are 
as follows: 


4 The stability of a system cannot be decided based upon the linearization 
when an eigenvalue of the linearized system is at the origin. 
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Fig. 7. Cathode inlet and the exit pressures, along with the stack pressure for 
oxygen excess ratio of 2. 


(1) Constant temperature is assumed throughout the stack, Ta = 
T, = Tin,a =de = 353 K. 

(2) Dry hydrogen is supplied to the anode inlet, Wy,0,in,a = 0. 

(3) A dead ended anode is considered, Wout,a = Wu>0,out,a = 
0, where all the fuel supplied is consumed and the hydrogen 
excess ratio equals 1. 

(4) The air flow rate is selected to provides the desired oxy- 
gen excess ratio” (approximately 2) required for an optimal 
operation of the stack [12]. 

(5) The air supplied to the cathode is completely humidified, 
Pin,c =]. 

(6) The cathode inlet and exit pressures are shown in Fig. 7 
along with the pressure inside the stack. The stack model 
parameters necessary for the simulations are provided in 
Appendix, Table 2. 

(7) Water flowing out of the cathode is only in vapor phase, 
Wh20,1,0ut,c =0, 


Due to the last assumption, the conditions encountered in this 
section belong to category (a) mentioned in Table 1. 

These operating conditions are chosen as an example to 
demonstrate the model prediction and investigate the predicted 
flooding and drying steady state patterns for various loads (cur- 
rent drawn). In the following sections, the stack hydration char- 
acteristic, that is, the degree of flooding or drying in the anode 
and cathode electrodes, is quantified using the water activities 
in the two electrodes. Note that, the water activity, aw ,.), is used 
because it is proportional to the mass of water in the electrode 
(see Eq. (7)). When aw,() < 1, the water vapor pressure in the 
electrode is below saturation pressure (sub-saturated condition). 
Conversely, when dy.) > 1, the water vapor pressure in the 
electrode has reached its saturation value and liquid water is 
present (saturated condition). The predicted steady state values 


5 Oxygen excess ratio is the ratio of rate of oxygen supplied to the cathode to 
the rate at which oxygen consumed in the electrochemical process. 
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Table 2 

Values of parameters used for simulations 

Symbol Parameter Value 

tm Membrane thickness 1.275 x 107? m 

Mm, dry Molecular weight of dry membrane 1.1 kg mole7! 

Pm,dry Density of dry membrane 2 x 10° kg m~? 

Afe Active area of stack 0.028 m? 

Neells Number of cells in stack 381 

Tet Stack temperature 353K 

F Faraday’s constant 96485 Coulombs 

Va Anode volume 0.005 m? 

Vv Cathode volume 0.01 m3 

M50 Molecular weight of water 18.02 x 10-3 kg mole~! 
Mu, Molecular weight of hydrogen 2.016 x 107? kg mole! 
Mono Molecular weight of dry air 29.00 x 107? kg mole~! 
ap Scaling on diffusion coefficient 5 

kine Cathode inlet flow coefficient 0.36 x 1075 kg s7! Paw! 
Kout.c Cathode outlet flow coefficient 0.22 x 1075 kgs7! Pa! 


of anode and cathode water activities, a&w,c and aw,a, and their 
corresponding water contents, A, and Àa, respectively, are shown 
in Fig. 8 for different stack currents, Ist. The absolute value of 
stack current is indicated by the lower horizontal axis while the 
upper axis is used to denote the corresponding current density, 
İst = Ist/Afc. 

At all current densities, the electro-osmotic drag creates a wa- 
ter flow from anode to the cathode, which depletes the mass water 
in the anode while increasing the water in cathode. This creates 
a gradient between the water activities on the anode and cathode 
sides of the membrane. The gradient causes a back-diffusion 
from the cathode to the anode, whose magnitude depends upon 
Àc — Ag, the difference in the water contents on both sides of 
the membrane (see Eq. (3)). For the dead ended anode with dry 
anode inlet flow (i.e., Wx,0,in,a = WH>0,out,a = 0) considered 
here, anode side water equilibrium, dmy,0,a/dt = 0, is attained 
when A. — A, is such that the back diffusion balances the electro- 
osmotic drag, that is, Wmemb = 0. Therefore, for anode water to 
be at equilibrium, the conditions A, < àc, and correspondingly, 
dw.a < Gw,c, must be satisfied as seen in Fig. 8. 
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At the low current densities (ist < 0.2 Acm~ for the fuel 
cell model considered here), the cathode inlet flow humidifi- 
cation and water generated are inadequate to keep the cathode 
completely humidified, and hence, at steady state the cathode 
is sub-saturated, that is, ay,¢ < 1. The mass of water inside the 
anode also reaches steady state such that aw,a < aw,e < 1. This 
condition belongs to case (i)-(a), and as discussed in Section 3, 
it corresponds to water equilibria with stable dynamic behavior. 
The physical mechanism responsible for an increase in cathode 
humidity with stack current is as follows. The rate of water gen- 
eration increases with stack current, and is accompanied by an 
increase in mass of water inside the cathode, my,0,c, and the 
water activity, dy¢, for given oxygen excess ratio. The increase 
in MH,O,c Causes a higher mass fraction of vapor in the gas leav- 
ing the cathode and allows sufficient increase in Wy,0,out,c to 
achieve water balance in the cathode, dmy,0,</dt = 0. 

Once aw, = 1, that is, mH,0,c = mibo, additional in- 
crease in mass of water in the cathode does not affect Wp,o,out,c 
for a given flow rate of dry gas, due to the assumption that 
only water vapor is carried out of the electrode. Above ist = 
0.2 Acm7?, the water removed from a saturated cathode is 
less than the total water entering the electrode through wa- 
ter generation and humidification of cathode inlet flow. There- 
fore, dmy,0,c/dt > 0, and excess water accumulates as liquid 
in the cathode while a part of water is transferred to the anode 
due to back diffusion through the membrane. The condition at 
ist = 0.2 A cm~? corresponds to an unstable equilibrium of case 
(iii)-(a), and my,0,c and my,0, increase with time due to an im- 
balance between the water supply and removal from the stack. 
This is shown in Fig. 8(a) with the vertical line at Zst = 60A 
(ist = 0.2 Acm7?). 

If the fuel cell operates long enough at ist = 0.2 A cm7?, the 
mass of water in the cathode exceeds a value such that, aw,c > 3, 
in which case the cathode water content saturates at 16.8 as 
shown in Fig. 1. This condition corresponds to case (v)-(a), in 
which the membrane water transport is independent of the mass 
of water in cathode, and only depends on the anode water con- 
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Fig. 8. The predicted anode and cathode water activities and water contents at steady state without water management. The cathode conditions are shown by the 
solid black lines. The anode side conditions depend upon the membrane characteristics and are shown for two different diffusion coefficients. 
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tent. This dependence provides the stabilizing factor for aya. 
The mass of water in the anode equilibrates such that the wa- 
ter diffusion through the membrane from cathode to anode is 
balanced by the electro-osmotic drag. As shown in Fig. 8, the 
discontinuity in aw,a, is a direct consequence of the jump in 
cathode water content (from 14 to 16.8) at the current density 
where the stack begins to flood. 

The increase in stack current causes a higher electro-osmotic 
drag which tends to reduce ay and increase dy c. The resulting 
increase in the difference between the anode and cathode side 
membrane water contents, A, — Aa, provides sufficient diffusion 
to balance the electro-osmotic drag, so that Wmemb = 0. Note 
here that, below ist = 0.2Acm~?, although A, — Aq increases 
with stack current, the increase in A, leads to an increase in 
Aa, and hence aw,a, with current. Above ist = 0.2 A cm, Ac is 
saturated at 16.8, and hence, increase in A, — A, is caused by the 
decrease in aw,a, and this decrease in aw,a continues for higher 
current densities (ig, > 0.2 Acm™~?). 

At sufficiently high stack current densities, the membrane 
water transport is inadequate to keep the anode completely hu- 
midified, and hence, aw,a < 1, leading to dehydration of the 
anode side of the membrane. The slope of the aw,a versus Ist, 
as well as the current densities where the anode starts to dry, 
depend upon the membrane water transport properties, namely, 
the coefficients for diffusion and the electro-osmotic drag. Ex- 
perimental data of the current densities at which we observe 
transitions from sub-saturated to saturated conditions [25], and 
back from saturated to sub-saturated conditions in the anode can 
be used to tune the membrane model and especially the coeffi- 
cients @p and ng. For the rest of the paper, ap = 5 is used such 
that the anode returns to sub-saturated conditions around 150A 
(ist = 0.54 Acm7?). 

The membrane model developed in this paper allows the pre- 
diction of the onset of anode flooding due to cathode flooding at 
low current densities and the anode drying at high current densi- 
ties [10]. These phenomena have been observed experimentally 
using various visualization techniques [25,2]. The flooding and 
water accumulation inside the stack further can be prevented by 
controlling flows entering and/or leaving the stack, as discussed 
in subsequent section. 


5. Stack water management 


The stack hydration characteristics predicted by the two- 
volume model are discussed in this section for two different 
water management strategies. One is using only cathode inlet 
humidity to achieve water management requirements, while the 
other is using maximum anode water removal plus appropriate 
cathode inlet humidification. 

The purpose of these water management strategies is to avoid 
water accumulation while maintaining saturated conditions in 
the anode and cathode electrodes, that is, dw.c, dw,a => 1. This 
constraint avoids sub-saturated conditions on the sides of the 
membrane that can decrease the stack performance [17]. Mean- 
while, least amount of liquid water in the anode and cathode is 
desired, since the presence of excessive liquid water inside the 
stack may block the reactant supply to the electrodes and affect 


the stack voltage. Considering that both sides of the membrane 
need to be 100% humidified, and aw,a < aw,c, this is achieved 
with ay, = 1. 

Without active water management, there is no guarantee for 
water balance and the water in the anode and cathode electrodes 
could deplete or accumulate, such as the case shown in Fig. 3. 
A control volume analysis of the fuel cell stack renders the rate 
of water accumulation inside the stack, Wy,0, stack, a8 


W450, stack = We 0,gen + WHO, inc — WH20,out,c — Wrem,a; 
(48) 


where for the ease of notation we define the net water removed 
from the anode side, Wrem,a, aS 


Wrem,a — WH) 0,out,a = WO, in,a- (49) 


Thus, the water management objective is to prevent water ac- 
cumulation inside the stack by maintaining Wy,0, stack = 0 and 
to regulate aw a = 1. These water management requirements 
should be met using appropriate cathode inlet humidification 
along with water removal from anode side, subject to the con- 
straint that the air flow is controlled to meet the oxygen excess 
ratio. 


5.1. Water management using variable cathode inlet flow 


humidity 


The first strategy uses only cathode inlet humidity to control 
the rate of water entering the cathode such that the water manage- 
ment requirements are met without anode water removal, (that 
is, Wrem.a = 0). The desired cathode inlet water flow, Wit in, és 
to achieve Wy,0, stack = 0 is obtained by rearranging Eq. (48) as 


des — sat 
W550. in,c = ™H>0,out,c Who, gen; (50) 


where Wy,0,out,c is replaced by Wino, out,c> the water removed 
out of the cathode under saturated conditions (aw c > 1). Note 
here that, we impose condition (50) to find the requirements for 
achieving a fully humidified cathode. 

In Fig. 9(a), WERO, out,c is shown by the dashed line while 
W4,0,gen is shown by the dash-dot line for different stack cur- 
rents. The corresponding Wie in, c 1s given by the solid line in 
Fig. 9(b). In the figure, the dashed line represents WEBO in, x the 
maximum cathode inlet water flow in vapor phase that can be at- 
tained for a given dry gas flow rate (oxygen excess ratio). Since 


Wio, in, x corresponds to 100% humidified cathode inlet flow, 
the desired cathode inlet humidity, ø% , is then given by 


in,c? 


des 

des ss H20,in,c 

in,c min wsat 2 1 $ (S 1) 
H20, in,c 


which is shown in Fig. 10. 

Clearly at low current densities (ist < 0.2 Acm72), the de- 
sired water supply to keep the stack humidified cannot be met 
even with 100% humidification. Due to inadequate water supply 
to sustain dy,¢ = 1 at steady state, the cathode and anode elec- 
trodes will be sub-saturated as seen in Fig. 8. In this case, the 
amount of water vapor supplied to the cathode can be increased 


A.Y. Karnik et al. / Journal of Power Sources 164 (2007) 590-605 601 


Current Density (A/cm?) 
014 021 029 036 043 05 057 064 
T T T T 


Water Flow rate (1 ò? kg/s) 
oa 
‘ 
s 


40 60 80 100 120 140 160 180 
(a) Current Drawn (A) 


Current Density (A/cm?) 
014 021 029 036 048 05 087 064 
T T T T T 


Cathode inlet water flow (10° kg/s) 


2 L 4 1 4 4 L 
40 60 80 100 120 140 160 180 


(b) Current Drawn (A) 


Fig. 9. Flows that influence the water management requirements, and hence, the desired cathode water inlet flow rates required for stack water balance at various 


current densities. 


by modifying the oxygen excess ratio, which is beyond the scope 
of this paper. As the current density increases, the increase in 
water generated through the electrochemical reaction reduces 
the humidity requirement for the cathode inlet flow, as shown in 
Fig. 10. 

Fig. 11(a) shows the steady state cathode and anode water 
activities obtained by using this desired cathode inlet humidi- 
fication, and regulating anode water activity at 1. The cathode 
inlet humidity is set using a feedforward map from the current 
as shown in Fig. 11(b) to provide the water balance. The feed- 
forward map, however, will not guarantee that the desired water 
activities are attained, since the desired equilibrium conditions 
belonging to case (iii)-(a) are unstable. The desired anode water 
activity is obtained by using a feedback control scheme, such as 
the one shown in Fig. 11(b). 
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Fig. 10. Desired cathode inlet flow humidification to maintain Wp, ,O,stack = 0 
with no water removal from anode side. 


As the current density increases, the higher difference re- 
quired at steady state between A, and àa is achieved with an 
increase in aw,c. Since increasing dy,- above 3 does not in- 
crease the diffusion, therefore, beyond ist = 0.54Acm~2 the 
anode cannot be kept completely humidified due to the limita- 
tions imposed by the membrane water transport characteristics, 
and maintaining a water balance will result in aw,a < 1. In this 
case, the maximum anode humidity can be achieved by main- 
taining dy.¢ = 3. 

Although controlling the cathode inlet humidity could clearly 
achieve the stated requirements between ist = 0.2 Acm~? and 
ist = 0.54Acm7?, lowering cathode inlet humidity is not al- 
ways desirable. First, controlling the inlet humidity fast and ac- 
curately is limited by the slow temperature dynamics of the hu- 
midifier [26]. Second, reducing the cathode inlet humidity can 
dry the cells closer to the air inlet manifold and thus increase 
cell to cell variations. 

Alternatively, anode water removal can be employed in coor- 
dination with cathode inlet humidification to achieve the desired 
water management requirements, as investigated in the next sec- 
tion. 


5.2. Water management using anode water removal with 
variable cathode inlet humidity 


The water management requirement delineated at the begin- 
ning of this section can also be achieved using water removal 
from the anode [15]. If anode water removal alone is inade- 
quate to satisfy the water management requirements, only then 
the cathode inlet flow humidity is reduced below 100% to com- 
plement the anode water removal. The water removal can be 
obtained using excess hydrogen fuel, which is typically desired 
in order to achieve adequate fuel supply at the downstream end 
of the anode electrode. If this excess fuel is being recirculated 
[27,14], then a water removal device can be used to remove the 
water from the anode exit stream before the recirculated flow is 
introduced in the anode inlet side of the fuel cell system. 
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Fig. 11. Anode and cathode steady state response using desired cathode inlet humidification at steady state to maintain Wy, 0, stack = 0 is shown in (a). In (b), the 


control architecture used to regulate the anode water activity to 1 is shown. 


The required anode water removal leading to WH,0, stack = 0, 
can be expressed as 


e sat 
Wretia = WhO, gen + WH 0, in,e — WEBO, out,o> (52) 


by rearranging Eq. (48), where Wino, out,c fepresents water re- 
moved from the cathode when the cathode is saturated, a condi- 
tion that is satisfied for ig. > 0.2 A cm7? (see Fig. 8). 

Fig. 12 shows the anode water removal required for differ- 
ent cathode inlet flow humidities. The maximum anode wa- 
ter removal is required with saturated cathode inlet flow, that 
is, when Wy,0,in,c = Wiino.in,c: Despite saturated cathode in- 
let flow, a negative anode water removal is required when 
ist < 0.2 Acm™, since the water generated and supplied are 
not enough to maintain saturated condition in the cathode, as 
shown in Fig. 8. Instead of anode water removal, the anode in- 
let flow should be humidified when Waem,a < 0, to maintain a 
saturated stack. 
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Fig. 12. Anode water removal required to maintain Wy, 0, stack = 0, for different 
humidities of cathode inlet flow. For conditions with Wrem,a < 0, additional 
humidification is required to keep the both sides of the membrane saturated. 


We now establish the feasibility of the anode water removal, 
that is, for which stack currents is it possible to remove all this 
water from the anode. The anode water removal capacity is lim- 
ited by the amount of water transported through the membrane. 
For aw,a > 1, the maximum water transported from the cathode 
to the anode is achieved with the maximum difference in water 
contents on the cathode and anode sides. This maximum wa- 
ter transport from cathode to anode, aera is achieved when 
aw,a = 1 with Ag = 14, and Ac = Amax, for which aw,c > 3. 

A thinner or a more diffusive membrane can provide a greater 

memb: 10 Fig. 13(a), Wine, corresponding to p = 5 is shown 
by the dashed line. The decrease in W)"*,, with increase in cur- 
rent drawn is due to the increase of water transferred from the 
anode to the cathode as the result of increased electro-osmotic 
drag. In the same figure, Wrethas the anode water removal re- 
quired to prevent the fuel cell stack flooding with 100% humid- 
ified cathode inlet flow, is shown by the dash-dot line in Fig. 
13(a). 

Above ig = 0.43 A cm7?, the rate at which the water accu- 
mulates in the fuel cell stack is higher than the maximum rate 
at which water is transferred to the anode side, thereby limiting 
the authority of anode water removal. Based upon Fig. 13(a), the 
desired anode water removal, Ws ,, to prevent anode flooding 
or drying is given by the minimum of the two values between 


the dash-dot and the dashed lines, that is, 


des _ JV: re max 
Wrem,a = min 1 Wiech ae W ahb] * (53) 
This W{°s , meets the anode water management requirement, 


but will not be able to meet the water removal requirement for 
the stack for ist > 0.43 A cm7?. In this case, the desired cathode 
inlet flow humidity, gees necessary to achieve Wy,0, stack = 0, 


: h in,c? 
is given by 
des sat 
des _ Wrem,a F Wibo, out,c 7 Wi, gen 1 4 
Pinc = min wsat , (54) 
H20,in,c 


The desired cathode inlet flow humidity decreases with increase 
in ig. beyond 0.43 A cm7?, as shown in Fig. 13(b). For com- 
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Fig. 13. Water management using anode water removal and cathode inlet flow humidification. In (a), the desired anode water removal is the minimum of dash-dot 


wmax 


line, representing Wye, and dashed line, representing enb 
parison, the desired cathode inlet humidity for the system with- 
out anode water removal is also included in Fig. 13(b) as the 
dash-dot line. The use of anode water removal allows cathode 
inlet air to be fully humidified in a wider current density range 
(0-0.43 Acm~), compared to the system without anode water 
removal (0-0.2 Acm~”). 

Beyond ist = 0.54 Acm~2, the maximum membrane water 
transport from cathode to anode is negative when the anode is 


100% humidified, hence, even with Wee = Othe anode will be 
sub-saturated (awa < 1). With wees a = 0, water management 


strategy similar to that discussed in Section 5.1 should be used, in 
which, only cathode inlet flow humidification is used to mitigate 
stack flooding. As seen in Fig. 13(b), the desired cathode inlet 
humidity is identical for both cases when ist > 0.54 A cm”. The 
humidification of anode inlet flow is then necessary to maintain 
awa = 1l. 

Fig. 14(a) shows the cathode water activity obtained by using 
desired anode water removal and desired cathode inlet humid- 
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ification. As discussed earlier the anode and cathode are sat- 
urated, and hence, have unstable water equilibria in the stack. 
For 0.2 A cm7? >i st > 0.43 A cm? the equilibrium conditions 
belong to case (iii)-(a) discussed in Section 3, the anode wa- 
ter removal is controlled using a feedback controller to main- 
tain the aw, = 1, as shown in Fig. 14. With 0.43 Acm—2> 
ist > 0.54 A cm, however, the desired equilibrium conditions 
belong to case (v)-(a) where the dynamics of anode and cath- 
ode are decoupled (see Fig. 6). Now, the mass of water in the 
cathode has unstable equilibria, and hence, the cathode water 
activity needs to be regulated when ist > 0.43 Acm™?. 
Compared to the strategy with only variable cathode inlet hu- 
midity, this strategy results in a higher steady state ay, between 
ist = 0.20 A cm7? and ig = 0.54 A cm7°, in order to support 
a higher diffusion through the membrane. Higher diffusion is 
required in this case to balance both electro-osmotic drag and 
anode water removal at steady state, as seen from Eq. (10). The 
advantage, however, is that the cathode inlet humidity can be 
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Fig. 14. Anode and cathode steady state response using desired anode water removal along with cathode inlet humidification at steady state to maintain Wy, 0, stack = 0 
is shown in (a). The control architecture used to regulate anode water activity to 1 for 0.2 A cm~? < ig < 0.43 Acm7? is shown in (b). 
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higher to avoid local drying at the cathode inlet. In addition, this 
system relies less on the cathode inlet humidity control, which 
often presents a challenge given the slow dynamics of the humid- 
ifier. It should be noted that even with anode water removal, the 
water management requirements can not be met for stack cur- 
rent densities below 0.2 A cm7? or above 0.54 Acm~2. In fact, 
this strategy results in the same anode and cathode activities as 
using cathode inlet humidification only. 


5.3. Summary of water management strategies 


The water management strategy discussed above maintains a 
saturated anode using maximum possible anode water removal 
and complementing it with appropriate cathode inlet humidity. 
This strategy is not unique. For example, the water management 
requirements can also be met by using a lower cathode inlet 
humidification accompanied by a lower anode water removal. 

The use of different combinations of anode water removal 
and cathode inlet humidity is summarized pictorially in Fig. 15, 
where each circle or square denotes a unit flow rate of water. The 
squares shown in the figure should be ignored in the absence of 
anode water removal. With only variable cathode inlet humidity, 
the sum of the circles shown in the first and second column 
(Wu30,gen and Wy, 0,in,c), Should equal the solid circles in the 
third column (W4,0,out,c) in order to achieve water balance. 
An excess of solid circles in the third column results in a sub- 
saturated cathode due to excess cathode water removal. 

The squares in the second column (Wy, 0, in,c) and fourth col- 
umn (Wrem,a) should be included when using anode water re- 
moval. In this case, the sum of circles and squares in the first two 
columns should equal the sum of circles and squares in the last 
two columns. Here, the additional cathode inlet humidification 
is balanced by providing appropriate anode water removal. 

Besides water balance, other factors may also constrain the 
desired anode water removal and cathode inlet humidity. For 
example, the excess hydrogen supply on the anode side, essen- 
tial for anode water removal, may be subject to constraints of 
specific recirculation device used. On the other hand, a conse- 
quence of combining anode water removal with higher cathode 
inlet humidification is the increase in cathode flooding when the 


Water Cathode inlet Cathode water Anode water 


Generated humidification removal removal 
+ + = A eee 
‘Anode & i 
Ho B O E 
Subsaturated 
eons: 
43 
5 (X) 
E + 8 &8 o 
O 
=+8 B W H 
[8] 
§ 
N 
+& > 8% u 
Sl Or ‘Anode ai | 
= Subsaturated 


Fig. 15. Summary of the different water management strategies required to avoid 
fuel cell stack flooding and specifically maintaining anode humidity at 1. 


anode humidity is regulated at 100%. These constraints need to 
be considered in specifying the requirements for anode water 
removal and cathode inlet humidity [16]. 


6. Conclusion 


A simple fuel cell stack hydration model is presented in this 
paper and is used for water management studies. The model pre- 
dicts anode flooding at low to medium loads once the cathode 
floods. The amount of anode flooding reduces with increase in 
current density and the anode exhibits drying at high current 
densities while the cathode is flooding. We show how the pre- 
dicted pattern is influenced by model coefficients such as the 
membrane water diffusion. This coefficient can be used to cali- 
brate the model and match experimentally observed patterns of 
electrode flooding as a function of the current drawn from the 
cell. 

The water management strategy necessary to maintain a hu- 
midified anode without the presence of any liquid in the anode 
is systematically analyzed. The feasibility of meeting the wa- 
ter management requirements with a combination of anode wa- 
ter removal and cathode inlet humidification is discussed for a 
high pressure fuel cell system without anode inlet flow humid- 
ification. The water removal requirement increases with stack 
current, however, the effectiveness of removing water from the 
anode side is limited by the water transport capability of the 
membrane which is restricted at higher current densities due to 
the large electro-osmotic drag. Therefore, at high current levels, 
reduction in cathode inlet humidification is necessary to meet 
the water management requirements. 

In addition to selecting the set-points for anode water removal 
and cathode inlet humidification, a feedback controller is nec- 
essary to regulate the system hydration at its desired levels. The 
need for a feedback controller is demonstrated by performing a 
stability analysis for the system under different operating condi- 
tions, and showing a lack of self stabilizing characteristic of the 
system at selected operating conditions. The stability analysis 
results depend upon the nature of interaction between the liquid 
water inside the electrodes and the water leaving the electrodes. 
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